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Abstract. In this paper we demonstrate for the first time that it is possible 
to solve numerically the Cauchy problem for the linearisation of the general 
conformal field equations near spacelike infinity, which is only well-defined in 
Friedrich's cylinder picture. We have restricted ourselves to the "core" of the 
equations: the spin-2 system, here, and as a background space-time we choose 
Minkowski space. We compute the numerical solutions for various classes of 
initial data, do convergence tests and also compare to exact solutions. We 
also choose initial data which intentionally violate the smoothness conditions 
and then check the analytical predictions about singularities. This paper is 
the first step in a long-term investigation of the use of conformal methods in 
numerical relativity. 



1. Introduction 

The subject of numerical relativity has made significant strides since 2005, when 
Pretorius [1] presented the first simulation of a binary black hole system, which 
extended over several orbits. Since then the numerical relativity community has 
developed more and more sophisticated codes to track the evolution of highly re- 
lativistic binaries of compact objects, increasingly including neutron stars. The 
codes are stable and accurate so that they can be used for production runs. In 
fact, the problem nowadays is not so much within an individual run but lies in 
the vastness of the parameter space (eight dimensional in the case of two spinning 
black holes). This calls for the development of sophisticated methods to exhaust 
this space [2]. 

However, there are still several problems with the numerical simulations, which 
have not been entirely solved. We mention here only the issue of the outer bound- 
ary and the feasibility of a global numerical simulation of an entire space-time. 
These are largely conceptual issues, which have been dealt with mostly by practical 
considerations. The background behind these issues is that all numerical simula- 
tions of binary systems are based on an idealisation: one assumes that the system 
under consideration does not interact with anything else, i.e., that it be an isolated 
system [3J 0]. 

Mathematically, this idealisation leads to the requirement that the space-time 
manifold describing the system be asymptotically flat. The main consequence of 
this is that the gravitational radiation, which is emitted by the system can be read 
off uniquely at null-infinity, commonly denoted by <f . This is a non-trivial result 
since the background independence of the theory (a consequence of the general 
coordinate invariance) prohibits the introduction of any absolute structure to which 
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the gravitational waves (the proverbial 'ripples' on space-time) can be referred to. 
Furthermore, the 'infinity' of an asymptotically flat space-time is unique. As a 
consequence gravitational radiation is a gauge dependent concept except at infinity. 
Note, that infinity in the space-time context refers to infinitely large distances as 
well as infinitely long times. So, the unique gravitational wave information from a 
radiating system can be detected only at infinitely distances from the system after 
infinitely long times. 

Obviously, this imposes a huge burden on the numerical simulations, which can- 
not incorporate infinite resources and, therefore, need to recourse to practical meas- 
ures. This is done by imposing an outer boundary, which is considered to be 'prac- 
tically' at infinity and the gravitational field is traced along the time-like hyper- 
surface swept out by the boundary during the simulation. The field is suitably 
interpreted as radiation data using appropriate 'extraction schemes'. 

The introduction of an artificial boundary into an infinite problem usually comes 
with several problems. Not only is it necessary to specify boundary conditions but 
these should be physically reasonable and, in addition, should lead to a mathemat- 
ically well-posed initial boundary value problem (IBVP) and a numerically stable 
evolution scheme. The issue of the introduction of a time- like artificial bound- 
ary is physically relevant, as was pointed out by Zenginoglu [5], since the decay 
rates of the tails are different on time-like boundaries compared to null-infinity. So 
one cannot simply replace one with the other. In any case, introducing an artificial 
boundary is always an additional approximation beyond the conceptual idealisation 
of an isolated system. 

Experience shows that for all practical purposes the current codes are good 
enough. However, from a puristic point of view the situation is still unsatisfactory. 
In this paper we will present the first step of a numerical approach, which will 
ultimately allow us to evolve space-times globally from asymptotically Euclidean 
initial data. Our method is based on the geometric description of asymptotically 
simple space-times developed by Penrose and their analytical treatment pioneered 
by Friedrich. 

In 1965, Penrose [HI El had realised that asymptotically flat (and, more generally, 
asymptotically simple) space-times can be characterised by the fact that they have 
a conformal boundary, i.e., they can be isometrically embedded into a larger Lorent- 
zian manifold after a conformal rcscaling of their metric such that the boundary 
of the embedded manifold is given by the zero-set of the conformal factor. This 
means that these boundary points must be regarded as being at infinite distance 
from any event in the original 'physical' space-time. 

Of course, such a geometric characterisation does not imply anything about the 
existence of general solutions of the field equations, which do in fact admit such 
a conformal boundary. This gap has been filled by Friedrich [51 [9] using what 
he called the conformal field equations (CFE), a slight generalisation of Einstein's 
field equations. One of his main results was the theorem on the well-posedness of 
the hyperboloidal initial value problem. This means that if appropriate data on a 
hyperboloidal initial hyper-surface are given then there exists a solution of the CFE 
at least for some time in the future of the initial hyper-surface. This result implies 
that the corresponding physical space-time admits a smooth conformal boundary, 
i.e., it is asymptotically flat. If the initial data are close to Minkowski data then 
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the solution exists for such a long time that the physical space-time is in fact future 
complete, i.e., it exists for all times. 

This semi-global existence result has been combined with a result by Corvino fTO] 
on the existence of asymptotically Euclidean initial data with an exact Schwarz- 
schild exterior by Chrusciel and Delay |TT] to obtain a global existence theorem 
for asymptotically flat space-times. Thus, Penrose's conjecture on the existence of 
general asymptotically flat space-times has finally been settled to the affirmative. 

The CFE operate on the larger, conformally related 'unphysical' manifold, within 
which the original 'physical' space-time occupies a finite set. This property makes 
it very tempting to use the CFE for numerical simulations because they capture the 
entire physical space-time without the need for an artificial boundary. The feasib- 
ility of this idea has been shown in several different scenarios [T2l HUI ITU [TBI [T51 [TT1 
1151 rn?l |2"U1 |2"T1 The main advantage of this approach is, of course, the fact 

that the global physical space-time is covered by one computational domain. It is 
worth mentioning that we are dealing here with a conformal rescaling of the metric, 
so that all the causal relations and the propagation properties of the (gravitational) 
waves are identical in the physical and the unphysical manifolds, see [24]. Further- 
more, since the conformal boundary is accessible from the unphysical space-time 
this implies that the unique radiation data produced by the system can be read off 
directly on the boundary without any additional extraction procedure. 

The main drawback of this approach (apart from the increased number of un- 
knowns) has been the fact that the conformal factor itself is a variable in the system 
and, therefore, evolves. Hence, its zero-set is not known beforehand (unless one em- 
ploys a particular shift gauge called ' J^-freezing' [H]). It has to be located within 
the computational domain at each time-step in order to determine the radiation 
data. 

Another issue has been the fact that initial data need to be specified on hyper- 
boloidal hyper-surfaces. This has the consequence that only prediction is possible 
but not refrodiction. Also within this hyperboloidal initial value problem it is not 
possible to study effects of gravitational radiation coming in from null-infinity, such 
as for example the scattering of gravitational waves. For such questions it is ne- 
cessary to include the region around space- like infinity, the avoidance of which has 
been the motivation for introducing hyperboloidal hyper-surfaces in the first place. 

In [25 Friedrich describes another set of conformal field equations, the gener- 
alised CFE (in contrast to the earlier 'metric' CFE), which can be used to study 
exactly the region near space-like infinity. These GCFE make full use of the con- 
formal geometry not only on the level of the metric but also on the level of the 
connection. This has the important consequence that it is now possible to fix the 
conformal factor beforehand as a function of the coordinates so that the coordinate 
location of J? is known at every instant. Friedrich also shows that one can 'blow 
up' space-like infinity so that it is represented by a cylinder connecting past and 
future null-infinity. 

The main advantage of this new system of conformal field equations is the fact 
that one can employ a gauge, such that almost all the equations take the form 
of simple advection equations along the time- like coordinate vector, i.e., ODEs. 
The only equations not of this type provide the evolution of the rescaled Weyl 
curvature, what Penrose called the gravitational field. This system is obtained 
from the Bianchi identities and it is easily shown to be symmetric hyperbolic. 
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The cylinder is a total characteristic for the GCFE, i.e., the equations reduce to 
an intrinsic system of evolution equations on the cylinder, emphasising the sole 
purpose of the cylinder, namely to provide a connection between the very late 
ingoing radiation and the very early outgoing radiation. 

These properties of the GCFE call out for an attempt to use them for numer- 
ical simulations. The fixed conformal factor, the ensuing fixed location of in 
the computational domain, the fact that one can cover a neighbourhood of space- 
like infinity including and beyond and and the advantageous form of the 
evolution system are all too tempting. 

However, all this comes at a price. The system of equations intrinsic to the 
cylinder degenerates on the locations, where the cylinder meets '. This poses 
serious analytical problems because generic solutions will develop singularities there, 
which are expected to travel along null-infinity, thereby spoiling its smoothness and 
hence the concept of gravitational radiation. In addition, the numerical simulations 
might be affected by the degeneracies in the evolution equations. 

In order to get a feeling for the possible effects we start out in this paper with a 
simplified situation. We study the gravitational perturbations on Minkowski space- 
time near space- like infinity in the gauge, which exhibits the cylinder described 
above. This amounts to solving the linear spin-2 zero-rest-mass field equation 
on a particular conformally flat space-time. However, as pointed out above, this 
system already captures the essential properties of the full non-linear system. The 
equations show the exact same behaviour as the fully non-linear GCFE in the sense 
that they reduce to an intrinsic set of evolution equations on the cylinder and they 
degenerate on the intersection with null-infinity. Our goal is to find out how the 
numerics reacts to the above mentioned degeneracies and how one can control the 
solutions at these locations. 

The paper is organised as follows. In sec. [2] we present the system under study. 
In sec. [3] we describe some analytical properties of the system which are relevant 
for the solution. Sect. [4] contains the description of the numerical implementation 
and presents our results. We close the paper with a discussion in sect. [5] 

2. The spin-2 system near space-like infinity 

In this section we derive the spin-2 zero-rest-mass equation on Minkowski space- 
time near space-like infinity. This section follows largely the exposition in the work 
by Friedrich [2"51|2"5], 

2.1. Minkowski space-time near spatial infinity. The metric of Minkowski 
space-time M in Cartesian coordinates X a is given by 

(1) g = r lab dX a dX b = dT 2 - dR 2 - R 2 dcu 2 

where, rj ab = diag(l, -1, -1, -1), T = X°, R 2 = {X 1 ) 2 + (X 2 ) 2 + (X 3 ) 2 and dco 2 
is the unit sphere metric. The exterior M_ of the null-cone of the origin is given 
by X a X a < 0. In order to focus on space-like infinity we now perform a reflection 
at the origin on M_ by defining 
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This puts the metric into the form 

/0 s - r] ab dx a dx b 

( 2 ) 9 



(x ■ x) 2 

This metric is singular whenever x ■ x — 0, i.e., on the null-cone at infinity. So we 
define a conformally related metric 

(3) g' = fl 2 g = r] ab dx a dx b , Q = -(x-x) 

which extends smoothly to the null-cone of infinity. Note, that space-like infinity 
with respect to Af_ is represented in this metric as the point x a — 0. In order to 
exhibit the cylindrical structure referred to above Friedrich [25 performs a further 
rescaling of the metric using a function n(r) = r/i(r) where r 2 = (x 1 ) 2 + (a; 2 ) 2 + (a; 3 ) 2 
and /i is a smooth function with fi(0) = 1. The rescaled metric is defined by 

(4) g = ±g> = ~S = ~ ( d (-°) 2 - dr * - r 2 ^ 2 ) . 

The final step is the introduction of a new time-coordinate t by x° = n(r)t. This 
gives the final form of the metric 

(5) g = -2 (n 2 dt 2 + 2tKK'dtdr - (1 - tV 2 )dr 2 - r 2 duj 2 ) . 
Note, that this metric is spherically symmetric. 

2.2. The cylinder at space-like infinity. We write the metric Q in the form 

, 2 „ ™' , dr (1-tV 2 ) fdr\ 2 1 , 2 

6 g = dt 2 + 2t — dt ^ , -[— - ^du 2 . 

k r fi z \ r J ^ 

This shows that the surfaces (i, r) = const have a non- vanishing area for each value 
of t and all r > 0. Since /i(0) — 1 this also holds for r = by continuity. In the 
coordinates (t,r) the set r = is a cylinder with topology K x § 2 . 

The conformal factor used between the Minkowski metric and the metric §5§ 

is 

fl r 2 — t 2 K 2 r , 9 9x 

e = - = — r — = -(i - tV). 



K 



Let Af_ = {6 > 0} denote the part of the conformal space-time, which corres- 
ponds to Minkowski space-time M . The boundary of Af_ consists of three parts: 
we have J !± = {1 =F t/i = 0} and we denote by / the set {—1 < t < 1, r = 0}. Fur- 
thermore, let I ± = {t = ±l,r = 0} denote the regions, where I touches J^ ± , and 
7° = {< = r = 0} the intersection of the hyper-surface t = with 7. In the following 
we will use mainly two choices for the function fi, namely (i = 1 or fi — 1/(1 + r). 
For the first choice, J is located at t = ±1, while the second choice puts it onto 
t = ±(1 + r), see Fig. [TJ 

In order to get some insight into the structure near space-like infinity we intro- 
duce double null coordinates 

(7) u = nt — 7', v = nt + r, 

which puts the metric into the form 

1 A A 1 A 2 

g = —audi; -duj . 
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Infinity for n(r) =1 



Infinity for /u(r) =1/(1 +r) 



— - u = const 

— v = const 

— infinity 




— - u = const 

— v = const 

— infinity 




FIGURE 1. The neighbourhood of space-like infinity in two differ- 
ent (t, 7')-coordinate representations corresponding to /i(r) = 1 and 
M (r) = l/(l + r) 



Null-infinity is characterised by the vanishing of one of the null coordinates, the 
other one being non-zero (u ~ 0, v ^ on J !+ and u ^ 0, v = on J^~). Both 
coordinates vanish on the cylinder /. In Fig. [T] we also display the lines of constant 
u and v. Since these are coordinates adapted to the conformal structure of M_ 
we can see clearly, that the cylinder is 'invisible' from the point of view of the 
conformal structure. This structure is exactly the same as the representation of 
space-like infinity in the coordinates x a on M_, it appears like a point. However, 
the differentiable structures defined by the (it, v) and the (t, r) coordinates are 
completely different near the boundary r — 0. 

2.3. The spin-2 equation. We use the spinorial form for the spin-2 zero-rest-mass 
field <j}ABCD — 4>(abcd)> see Penrose and Rindler |27| . the conventions of which 
(the GHP-formalism) we adopt throughout this work. The field equation is 



(8) 



•^ABCD = 0, 



where V aa' is the spinorial connection defined by a general 3 + 1 metric g. Using 
the compacted spin-coefficient formalism |28l 2? j and defining the spinor field com- 
ponents cf>k in the usual way (0o — ^oooo, <Pi = 0oooii- ■ ■ ) we obtain the equation^] 



In the following equation only, the symbol k stands for the spin coefficient from the GHP 
formalism and is not to be confused with the function k introduced before (for the metric g in 
Eq. (JgJ ) . Since, as will see shortly, the spin coefficient re vanishes identically under our assumptions, 
we continue to use the symbol re for the function in Eq. |6J. Moreover, it should always be clear 
from the context whether a prime corresponds to the operator from the GHP formalism or to a 
radial derivative as before. 
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for k = 1 : 4 

t>0fc - SVfc-i = -kr'(f>k-i - (4 - k)K(/) k+1 + (k- l)a'4> k -2 + (5 - fc)p0 fe , 

pVfc-i - = -(5 - fc)T0/c - (fc - l)«V*-2 + (4 - k)a<t>k + x + kp'4>k-i- 

The metric g in Eq. ([6| is spherically symmetric. Using a null-tetrad, which 
is adapted to this symmetry, we know immediately that all spin-coefficients with 
non-zero spin-weight vanish. The only non-vanishing spin-coefficients are 

, 1 , 

We choose the real null-vectors of the null-tetrad as 

(io) rd a = -J= ((i - tK 1 ) d t + Kd r ) , n a d a = -^=((i + tK')d t -Kd r ), 

and complement them with a complex space-like null- vector ni a d a , which is tangent 
to the spheres of symmetry; one particular choice is 

V2 V sm6> 

where (#,</>) are the standard polar coordinates on § 2 . Then the relationship 
between the GHP-operators p and p' and the directional derivatives D and D' 
(see [57]) are mediated by the unweighted spin-coefficients e and 7, which are 

1 , 

£ = 7 = pK , 

2y/2 

for the metric g. Furthermore, we redefine the (3 and <3' operators so that they refer 
to the unit sphere. This is achieved by the replacement 

(11) 

these new operators act on a function 77, which is defined on the unit sphere and 
which has spin- weight s, by 

9?/ = m a d a r] — s cot 9 7/, 9'?/ = fh a d a r\ — s cot 8 rj. 
With these specialisations the eqs. (|9| become 

(1 - tK')dt(/)k + nd r (j) k - (3k' - (5 - fc)/i)0fe = /jB'<j>k-i, fc = l:4, 
(1 + tn')d t (t)k - nd r 4>k + (3k' + (k + l)fJ,)4>k = pd(f>k+i, k — : 3. 

Finally, we use the spherical symmetry of the background metric to expand the 
spinor component <pk into spin-weighted spherical harmonics s Yi m with s — 2 — k. 
\s\ < I and \m\ < I, i.e., we write for any (t,r, e) elx M+ x S 2 



00 1 



(t,r,e)=^ E ^k m (t,r) s Y lm (e) 

l = \s\ m= — l 



Using the properties of the unit-sphere 8 and 3' operators acting on spin-weighted 
spherical harmonics 

3sY lm = + s(s + l) s+1 Y lm , B' s Y lm = y/l(l + 1) - s(s - 1) s -iY lm , 
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we can decouple the equations and obtain a system of PDEs in 1 + 1 dimensions 
for any admissible pair (l,m) 

(1 - tn')d t( j> l r + Kd r cj>[ m = {3k' - 4fi)cj>[ m + M a 2 0d m , 

(1 - tK')d t <f> l 2 m + Kd r <f> l 2 m = (3k' - 3/i)4 m + /J,a (t> l r, 

(1 - tn r )d t 4> 1 ^ + K d r (j> l r = (3^' - 2^)4™ + M«o02"\ 

(1 - tn'^r + Kd r cj> l 4 m = (3k' - M )^ m + fia 2 <p l 3 m , 



(13) 



and 



(14) 



(1 + t K ')d t (f> l 3 m - nd r (t) 1 ^ = -(3k' - 4/i)^ m - \ia*l>Ti 
(1 + tK')d t (t> l 2 m - Kd r (f> l 2 m = -(3k' - 3jj)(j) 2 m - Aiao^" 1 



y 3 . 

(1 + tK')d t (t>[ m - Kd r <j>\ m = -(3k' - 2fj,)<t>[ m - fj,a (p 2 m , 
(1 + tK^dt^ - Kd r <t> 1 ™ = -(3k' - fi)^ - fia 2 <p[ m . 

We have defined the two quantities ao = \Jl(l + 1) and a 2 = \/l(l + 1) — 2. We 
also abuse notation by omitting the superscript (I, m) at the spinor components. 
Thus, we have to keep in mind that we get a system like the above for every 
admissible pair (l,m). Note, however, that the for the cases I = and I = 1 the 
system looks different because a 2 = or = 0. These are advection equations 
along the null-directions defined by l a and n a , from which can deduce a system of 
evolution equations and constraint equations. 

3. Analytical properties 

Before we describe the numerical implementation we want to briefly discuss some 
of the relevant analytical properties of this system. 

3.1. Evolution and constraint systems. We first split the system into evolu- 
tion equations and constraints by taking appropriate linear combinations of the 
equations. This yields the evolution system 

(1 + tK')d t <fio - nd r (j>o = -(3k' - n)4> - fia^i, 

( 15 ) d t (j)2 = ^fJ-a <Pi - ^"003, 

dt4>3 = M03 + ^M a O02 - ^M a 204, 
(1 — tK')dt4>4 + Kd r (j)4 = (3k' — /i)</>4 + /i«203 

and the constraints 

Ci = —2Kd r (f>i + 6rfi'<j)i — 2tK'fi(j)i + a fi(l — tK')<p 2 + a 2 fi(l + tK r )(f> = 0, 

(16) C 2 ee -2Kd r <j) 2 + &r^'(t> 2 + a M(! - tn'Jfa + a /J,(l + tn')^ = 0, 

C3 ee —2Kd r (f>3 + 6rfi'<f>3 + 2tK r n4>z + «oA*(l + tK')<p 2 + a 2 fi(l — tK r )(f>4 = 0. 

Obviously, the evolution equations can be written in the form 

A°<9 t $ + A 1 ^ = B$ 

with symmetric matrices A and A 1 . Furthermore, the matrix A = diag(l + 
in' , 1, 1, 1, 1 — tK 1 ) is positive whenever \t\ < 1/\k'\. Thus, the evolution equations 
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are symmetric hyperbolic on the domain, where t is confined to these values. It 
is also straightforward to show that the constraints are propagated by the evol- 
ution equations, i.e.. they propagate solutions of the constraints into solutions of 
the constraints: taking time derivatives of the constraint quantities Ck, using the 
evolution equations and the constraints we obtain the 'subsidiary' system for the 
propagation of the constraint quantities 

d t C\ = -nCi - ~a ^C 2) 

(17) d t C 2 = ^aoniCt - C 3 ), 

d t C z = [J.C-3 + ^a [iC 2 - 

Since this is a homogeneous system of ODE, the claim follows immediately. 

3.2. Characteristics. It is instructive to study the behaviour of the characteristics 
for the evolution system. Clearly, we have three different characteristics, the integral 
curves of the vectors d t , l a d a oc (1 — tK r )d t + nd r and n a d a oc (1 + tn')d t — nd r . 
It is easy to see that the latter two are exactly the lines of constant u and v 
defined in Eq. [7] respectively, and hence, they are null. In Fig. [2] we show again 
the neighbourhood of / and ^ + for fi(r) = 1/(1 + r). This time we show the 
characteristics also in the unphysical part of the diagram. The non-shaded region 




Figure 2. The characteristics of the evolution equations in a 
neighbourhood of I + . The shaded region is the domain where the 
equations fail to be hyperbolic. The corresponding neighbourhood 
of I~ is obtained by reflection at the lower border of the diagram 
accompanied by an interchange of u and v. The thick (green) 
broken line indicates the boundary of the domain of hyperbolicity 
(shaded region). 



bounded partly by the thick broken line is the domain of hyperbolicity of the 
evolution equations. Notice that every neighbourhood of J + contains regions where 
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the hyperbolicity breaks down. Apart from the fact that one cannot hope to get 
existence and uniqueness of solutions beyond that region its presence also makes 
the numerical evolution challenging. For instance, setting up an initial boundary 
value problem with the left boundary at negative values of r does not make sense 
if the evolution is to reach up to / + . Even with the left boundary on I it is not 
possible to go beyond I + since the evolution hits the non-hyperbolicity region. 

3.3. The constraint equations. The constraint equations on t = are obtained 
from (16). Define ipk — H ^4>k and let Z > 2 so that ao and a 2 both do not vanish. 



Then 

2rd r ipi = a ip 2 + a 2 ipo, 

(18) 2rd r ip 2 = otQipi + a ip 3 , 

2rd r ip 3 = a o 02 + a 2 ip4- 

Solutions of this system can be obtained by specifying functions ipo(r) and 04 (r), 
giving initial conditions at r — and then solving the equations for the remaining 
components. However, for solutions with bounded derivatives at r — we cannot 
specify these data arbitrarily. Instead, for r = 0we have the equations 

= a o 2 (O) + a 2 o (O), 

(19) = a Vi(0) + a o 3 (O), 

= ao0 2 (O) + a 2 4 (O). 

For I > 2 this is an inhomogeneous linear system for the components 0i, 02, 03 
with inhomogeneity given by 0q an d 04- Since the homogeneous system has a non- 
trivial solution, we can obtain solutions of the inhomogeneous system only when 
0o (0) = 04 (0) holds. Then there is a 2-dimensional set of initial conditions for 
bounded solutions given by 

(20) 0o(O)=a, 0i(O) = 6, 2 (O) = - — a, 3 (O) = -b, 4 (O) = a, 

where a and b are arbitrary complex constants. If I = 1 then ipo an d 04 vanish, since 
they have spin-weight ±2, and this implies 02(0) = 0. For / = the spin- weighted 



functions -01 and -03 vanish and ao = 0. Then (19 1 docs not impose any condition 
on 02 (0). 

The procedure to obtain initial data just described is suggested on physical 
grounds. One specifies the 'wave-like' degrees of freedom represented by ipo an d 
04 and then computes from them the other components, which give rise to the 
energy-momentum. There is also another way to obtain initial data for the spin-2 
system: we specify 02 freely. When I > 1 we use (l8|2) to obtain ipi + -03 and we 
can specify the difference tpi — 03 freely. Then, when I > 2, we can use ([l8|l,3) to 
obtain 0o and -04. Proceeding in this way the regularity conditions at r = are 
automatically satisfied. It has the advantage that no differential equations must be 
solved. 



3.4. The total characteristic /. We go back to the full system of equations (12) 
and for the discussion in this section we specialise to the 'horizontal' conformal 
representation defined by /x(r) = 1. Evaluating the equations on /, where n{r) = 
yields a system of eight PDE, which contain no radial derivatives, i.e., they are 
intrinsic to the cylinder. Therefore, no information can flow into and out of this 
set — it is a total characteristic. The evolution of the spin-2 field on / is entirely fixed 
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by the initial data given at some initial time, say t = 0. The intrinsic equations 
'transport' the initial data up (and down) to the location where null-infinity meets 
the cylinder. 

It was proven in [25 that under the assumption of analyticity of the rescaled 
"unphysical" metric g near space-like infinity the components of the spin-2 zero- 
rest-mass field admit an expansion of the form 

oo ^ 

(21) Mt,r,e) = J2 -jXfc,p(*.e)rP, with s = 2 - k 

p=\s\ 

witr0 

p i 

(22) e) = E E <Pk? P {t)sY lm (e), 

l=\s\ m=-J 

where Xk,p(t> e ) = dr <f>k(t: r i e)| r =o- Note, that the number of Z-modes is bounded 
by the differentiation order p. 

The analytical behaviour of the solutions of these 'transport equations' near the 
cylinder is completely known, see [28, 29], so we will only briefly summarise the 
analysis here. After decomposing the components into spherical modes we obtain 
the expansions 

(23) r) = E -/k m p (t) r p , with s = 2 - k 

p=\s\ 

with $£(t) = #Vi m (*,0lr=0 and I < m. 

Let p > 2 since the cases p = 0, 1 are special. Taking radial derivatives of the 



equations ( |13[ ) and ( 14 1 and then restricting to r = we can derive coupled second 
order equations for the coefficient functions <f> ™ p (t) of the five components of the 
spin-2 field 

(24) (l-t 2 )^ m p +[2 S + 2(p-l)i] 0[™+[;G + l)-p(p-l)]4 m p -O, fc = 0:4, 

which is a Jacobi differential equation. From the discussion in [30] we find that 
this equation has regular polynomial solutions if and only if |s| < I < p. On the 



other hand, from (23 1 we see that I = p is an admissible term in the expansion so 
that the highest Z-terms in this expansion cannot be regular. In fact, one can show 
(see [30J ) that they contain logarithmic divergences at t = ±1. 
The solutions for I = p are of the form 

(25) 2 2 

+ 2^(1,1 + fl +p;2p + 2;—) 

for some constants A™ p and B™ p . The functions P2 p S ~ p S ~ P {t) are Jacobi polynomi- 
als, so they are regular, while the hypergeometric functions 2-F\(l, 1 + s + p; 2p + 
2; y^j) contain the logarithmic divergences. 

To simplify the notation we will drop the indices m and p from these coefficients 



because they are irrelevant to the following discussion. Inserting the solutions (25) 



2 Friedrich uses a slightly different class of functions for this expansion. Here, we use the 
closely related well-known spin-spherical harmonics. See app. [E]for a detailed discussion of the 
relationship between both systems of functions. 
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into the eight first order equations (131, (14) shows that the coefficients A k resp. 
B k are linearly related so that we can determine them all given only one of them, 
say A := A2 resp. B := B 2 : 



(26) 

and 
(27) 



A = A 4 



a a 2 



(p+l)(p + 2) 



A, A 1 = A 3 



B — B4 



a a 2 
Pip - 1) 



B, Bi = B 3 



ao 
P 



B. 



A, 



This implies that we can eliminate the logarithmic terms by putting B = 0. We 
can express this condition in terms of the initial conditions at t — as follows: let 



(28) 



= ^(0) = A k 2- 2 P + B k l k , fc = 0:4 



where / = 2^1(1, 1 + s + p; 2p + 2; 2). By using (26 1 and eliminating A between 



successive equations ( 28 1 we obtain the four relations 



(29) 



a 2 ^ + (p + 2)^= ( £ + 
«o t 2 + (p + l)^j = (f 2 + 
a ± 2 + (p + 1)^ = f/ 2 + 
a 2 0, + (p + 2)0 d - f /, + 



p + 


2 


p- 


1 


p + 


1 


p 




p + 


1 


p 




p + 


2 


p- 


1 



a 2 ao 



P 
ao 
P 

p 



'B 



P 



Choosing initial data such that the left hand side in any of these relations vanishes 



implies that all left hand sides in ( 29 ) vanish (since one can show that the terms in 



parentheses on the right hand sides never vanish) and that the solution for <^™(t) 
on the cylinder has no logarithmic divergences. 
Notice, that these relations hold for p > 2. 



Since one can show 



that the 



coefficients of r° and r 1 in (23 1 are necessarily free of logarithms the first appearance 
of logarithmic divergences can occur for p = 2. We will use this fact in sect. |4.5| 

By going back to the expansion of the spinor components one can interpret ( |29[ ) 
as differential equations on the initial hyper-surface t — 0, relating angular and 
radial derivatives of the spinor components evaluated at r — 0. It is a consequence 
of a truly remarkable result of Friedrich's [23 which shows that these relations are 
the components of the equations 



(30) 



D 



{A l B 



,Da 2 



D A 



u p -BcDEF) 







evaluated at the cylinder. Here, Dab denotes the intrinsic (space spinor) derivative 
operator on the initial hyper-surface induced by the metric g and Babcd is the 
linearised Cotton spinor, expressible in our context as 

Babcd = 2D(a E Q 4>bcd)e + ^ D^a E 4>bcd)e, 

where il is the conformal factor relating the flat Minkowski space metric to the 
metric 
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3.5. Boundaries and boundary conditions. In this paper we discuss the solu- 
tion of the spin-2 equations in a neighbourhood of the cylinder at r = for times 
t < 1. Having decomposed the equations into different modes due to the rotational 
symmetry we have a 1 + 1 dimensional system. Our spatial domain is the interval 
[0, 1] and we are interested in the time development from t = up to as close to 
t = 1 as possible. 

We need to discuss the boundaries at r = and at r = 1. At the 'outer' boundary 
at r = 1 there is one ingoing characteristic, corresponding to (f>o and an outgoing 
characteristic corresponding to </>4. So we need to specify 4>q at r = 1. All the other 



characteristics are tangent to the boundary. Since the subsidiary system ( 17 1 has 
no ingoing modes at r = 1 there are no further conditions on the specification of 
())q apart from corner conditions at (t,r) = (0,1). 

The other boundary, r — 0, is a total characteristic, i.e., there are no ingo- 
ing or outgoing modes. All characteristic directions are tangent to the boundary. 
Therefore, we do not need to do anything special there. 

Let us point out that it would be possible to put an artificial boundary at some 
location r < 0, where again there would be one ingoing and one outgoing mode, i.e., 
one boundary condition. However, as discussed in sect. |3.2| it would be impossible 
with this setup to reach t = 1 because of the breakdown of hyperbolicity in any 
neighbourhood near t = 1. 

Also, note that the 'outer' boundary at r — 1 (which is really the inner bound- 
ary from the point of view of the original Minkowski space picture) is an artificial 
boundary which exists only due to our conformal representation of the neighbour- 
hood of spatial infinity based on an involution of Minkowski space. If we had chosen 
a different representation such as the one based on the standard representation of 
Minkowski space in the Einstein static universe, then we could have eliminated this 
artificial boundary entirely and we could have studied the system globally on one 
computational domain. This will be discussed in more detail elsewhere. 

4. Numerical implementation and results 

4.1. Numerical algorithms. We use the method of lines to solve the equations 
numerically. Each space-like surface, {(t,r) : r E [0, 1]}, is represented as a grid with 
constant spacing Ar. The spatial discretisation is obtained by replacing the spatial 
derivatives with a summation by parts (SBP) operator, [3"T1 13"2]. SBP operators 
are finite difference approximations to derivatives that obey a discrete version of 
integration by parts, [31]. This ensures that the matrix, that represents the first 
order spatial derivative, is semi-bounded with respect to a suitable norm on a vector 
space, [3T] , and it allows for energy estimates to be calculated for our discretization 



of our system of equations, ( 15 ), by following the arguments used in the continuous 
case. 

We have used the minimum bandwidth, restricted full norm, fourth order (on 
the interior and third order on the boundary) SBP operator given in |311 page 65]. 
We have tried several other SBP operators from the papers |31l I3"3l I3~4] but have 
found that our choice gives the smallest absolute error over each grid point. 

The semi-discrete system thus obtained is then solved by the standard 4th order 
Runge-Kutta solver. The time-step At used in the RK4 method is computed either 
with a constant CFL factor C — At/Ar or by computing the maximal possible 
time-step from the information about the characteristic speeds available at each 



14 



F. BEYER, G. DOULIS, J. FRAUENDIENER, AND B. WHALE 



instant of time using the criterion that the numerical domain of dependence be 
larger than the analytical domain of dependence. This 'adaptive' time-step was 



particularly useful for runs of the cases with horizontal J? (n = 0, see sect. 2.1) 
because in this case one of the characteristic speeds diverge for t — >• 1. This has the 
consequence that any choice of constant CFL factor eventually leads to instability of 
the code near r = 1. The adaptive time-step allows us to approach t = 1 arbitrarily 
closely. Of course, the draw-back is that the time-step becomes arbitrarily small so 
that we cannot reach t = 1 exactly. 

The initial data are either specified by hand from exact solutions or by solving 
the constraint equations using one of the two methods discussed in sect. |3.3| When 
solving the ODE's we use the DOP853 method, a Runge-Kutta method of order 
8(5,3) due to Dormand and Prince as described in |35j . 

In line with the recommendation of [36 we implement the boundary conditions 
by using the SAT method (see |321 I3~4]). The equation enforcing the boundary 
condition, in semi-discrete notation, is 

4>o,i = , 1 , , ( «« (Q0o)j-(3k- - m) <f>o ti - Hia 2 (pi t i 
(31) L + tK i\ ' i = l:N 

- T , 8i, N (4> ,N - b r (t)), 

1 tKjy 

where <f>Q, <j>\, k, k' and fi are grid-functions, i.e., vectors of length N with 0o,i etc. 
being the i'th component. Q is an SBP operator represented as a constant N x N 
matrix acting on R and b r (t) is the free function given at the right boundary. The 
parameter r is suitably adjusted to enforce the boundary conditions in a stable way. 
For details we refer to |31l \M\ . 

We have implemented the numerical algorithms discussed above in a Python 
code, which relies on the numerical routines of the NumPy and SciPy packages. 

4.2. Code test with an exact solution. Here we will reproduce numerically the 
exact solution of the system (15l-(16l derived in appendix [A] This solution refers 



to the case with diagonal J' (n — \ see sect. 2.1|. For this we prescribe initial 



conditions of the form (set t = in ( 49 ) ) 



( 32 ) ^ = j^>^ = jYT^^ 2= (^^ 3 = I^ A4= aTW- 

We also impose, as discussed above, only a boundary condition for 4>q at r = 1: 

(33) MtA)= (2 I ^- 



Evolving the initial data ( 32 ) with a fixed time-step and computing the difference 
between the computed and the exact solution we obtain the convergence plots for 
the components <fio , <f>4 at time t = 1 as presented in Figures [3] and [4] 

Table [l] lists the log 2 of the absolute error in the L 2 -norm and the convergence 
rates at time t = 1 for 4>o,4>4. We find the expected overall 4th order convergence. 
However, note that this refers to the time t = 1, i.e., when the equation become 
singular at the left boundary. So we can reach t = 1 without loss of convergence. 
Running beyond t = 1. however, very quickly leads to instabilities and code crash, 
as it is to be expected from the loss of hyperbolicity of the equations near t = 1 on 
the left boundary. 
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Convergence of (p at time = 1 .0 



-35 

o 

.U -40 
60 

■2 -45 





50 gridpoints 
100 gridpoints 
:200 gridpoints 








400-gt.idpom.t 
80p gridpoint 
















y\ /-~~^ -\ 

\ 


















i-' 









































Figure 3. Convergence of 4>q to the exact solution at time t = 1. 



Convergence of <f) 4 at time = 1 .0 



a 



50 gridpoin 
100 gridpoin 
- ""^.200 gridpoin 
400 gridpoin 
800 gridpoin 


s — 

s — 

s 








s 
s 













































































-15 
-20 
-25 
-30 
-35 



Figure 4. Convergence of cf>4 to the exact solution at time t = 1. 



The convergence plots for (Fig. [3]) show that the inflow boundary at r = 1 
is very well behaved, similarly the total characteristic at r = shows no problems 
in the case of </>4, which propagates away from it towards increasing r. We do, 
however, see high frequency features at r = for 0o & n d at r = 1 for ^4. In 
both cases these are numerical artifacts. They indicate that the modes hitting the 
boundary from inside the grid are reflected due to numerical inaccuracies. We could 
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log 2 (||A|| 2 ) 


Rate 


log 2 (||A|| 2 ) 


Rate 


50 


-24.899275 




-11.322955 




100 


-29.057909 


4.1586 


-14.245253 


2.9223 


200 


-33.378371 


4.3205 


-17.326037 


3.0808 


400 


-37.781307 


4.4029 


-20.487582 


3.1615 


800 


-42.216813 


4.4355 


-23.689231 


3.2016 



Table 1. Absolute error A compared to the exact solution using 
the L 2 norm and convergence rates at time t = 1 for O , 04 . The 
calculation was done with a fixed time-step. The first two columns 
refer to 0o, the others to 04. 



eliminate them both by using small amount of artificial dissipation but have chosen 
not to. In the case of the 'outer' boundary they could be eliminated by using a 
different conformal representation. 

4.3. Code test with general initial data. We specify initial data using the 
algebraic method described in sect. |3.3| We specify 

2 (O, r) 



(4r/6 2 ) 16 (r-6) 16 0<r<6 
b < r < 1 



This gives a 'bump' centered at r — b/2 and we chose b = 4/5 for these runs. 
We also chose 0i(Q, r) = 03(0, r). Then all the spinor components are fixed by 
the constraint equations at t = 0, see sect. |3.3| They have compact support in 
the interval < r < 1 and in order to satisfy the corner conditions we choose the 
boundary condition 0o(£, 1) = 0. We run the tests with the two different conformal 
representations specified by n = and n — 1. Tab. [2] and tab. [3] show the results. 
In both cases we obtain the expected 4th order convergence in O . However, in 
the 'horizontal' representation we see a loss of convergence in the 04 component, 
which is due to the fact that the characteristics for 04 approach so that the 
characteristic speed becomes large. The decrease in the convergence rate might be 
due to fact that the code is leaving its stability region. 





log 2 (||A|| 2 ) 


Rate 


log 2 (||A|| 2 ) 


Rate 


50 


1.383744 




6.860341 




100 


-2.646956 


4.0307 


2.268950 


4.5914 


200 


-6.615024 


3.9681 


-2.951419 


5.2204 


400 


-10.692111 


4.0771 


-5.655223 


2.7038 



Table 2. Convergence rates of the absolute error A in the L 2 
norm at t = 0.96 for compactly supported initial data evolved in 
the 'horizontal' conformal representation with n — 0. The error 
is computed with respect to a higher resolution run with 800 grid 
points. The first two columns refer to 0o, the others to 04. 



Finally, in Tab. 4; we show the convergence rates for a higher mode (/ = 10) using 
compactly supported data in the 'diagonal' representation. 
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log 2 (||A|| 2 ) 


Rate 


log 2 (||A|| 2 ) 


Rate 


50 


1.447922 




2.228200 




100 


-2.508450 


3.9564 


-1.646880 


3.8751 


200 


-6.479996 


3.9715 


-5.645109 


3.9982 


400 


-10.558692 


4.0787 


-9.729219 


4.0841 



Table 3. Convergence rates of the absolute error A in the L 2 
norm at t — 1.0 for compactly supported initial data evolved in 
the 'diagonal' conformal representation with n = 1. The error is 
computed with respect to a higher resolution run with 800 grid 
points. The first two columns refer to fa, the others to fa. 





log 2 (||A|| 2 ) 


Rate 


log 2 (||A|| 2 ) 


Rate 


50 


-3.822906 




-1.422306 




100 


-7.773146 


3.9502 


-5.454578 


4.0323 


200 


-11.755335 


3.9822 


-9.437786 


3.9832 


400 


-15.836066 


4.0807 


-13.522334 


4.0845 



Table 4. Convergence rates of the absolute error A in the L 2 
norm at t = 1.0 for compactly supported initial data with I = 10 
evolved in the 'diagonal' conformal representation with n = 1. The 
error is computed with respect to a higher resolution run with 800 
grid points. The first two columns refer to fa. the others to fa. 



4.4. Constraint violation. Our next topic is the behaviour of the constraint 



quantities defined in ( 16 ). We use the same initial and boundary data as in sect. 4.3 
with the initial bump centered at r = 0.3. We take I — 2 and n — 1. At each time- 
step we evaluate the three constraint quantities. In Fig. [5] we show the L 2 norm of 
the three constraint quantities at selected instants of time in the interval [0,0.8]. 
The size of the constraint violation remains almost the same. This is clearly due 
to the fact that we are dealing with a linear system and also only for a short time 
interval. 

4.5. Transport equations. Now we are in a position to attempt a reproduction 



of Friedrich's analytical results on the cylinder for /i(r) = 1. Recall from sect. 3.4 
that the components of the spin-2 zero-rest-mass field 4>abcd admit an expansion 
of the form 

(34) fa(t, r,e)=J2 -Xk, P (t> e) r p , with s = 2 - k 



p=l«l 



wi 



th 



(35) X**(*,e) -EE A m P (t)sY lm (e), 

l=\s\ m=-J 



where \k,p — dr P ^ 4>k\r=o- The discussion in 4.5 showed that the time dependent 
coefficients (i) can contain logarithmic divergences at t = ±1 unless one restricts 
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Figure 5. Constraint violation for several values of t € [0,0.8]. 
We show the L 2 -norm of the three constraint quantities. 



their initial values at t = 0. The singular behaviour affects those coefficients with 
p — I and p > 2. 

In our numerical computation we have already decomposed the spinor compon- 
ents (£>k into their spherical harmonic pieces. So each function 4>k(t,r) in (13) 
and (14 1 already corresponds to a particular pair (l,m) and each of these modes 
has an expansion near r = of the following type 

oo p 

(36) Mt,r) = £ £ ^k m p (t)r p , with s = 2- k. 

P=\s\l=\s\ P - 

We choose the functions (/>o(0, r) and 04(0, r) at t — in such a way that 

$2 = lim^o(0,r)/rV0, 
and similarly for 4 (O,r). Specifically, we take 

(37) o (O,r) =8r 2 (r-l) 36 , 4 (O,r) = -8r 2 (r - l) 36 . 

With these functions we solve the constraints at t = and then we evolve the system 
for / = 2 and arbitrary m using the conformal representation with n = 0. We use 
the boundary condition <po(t, 1) = at the 'outer' boundary. The computation was 
done with the adaptive time-step in order to come as close to the singular time 
t = 1 as possible. It was stopped at t = 0.9999. 

On the other hand, following the procedure outlined in sect. [43] we can compute 
the exact form of the singular coefficients along the cylinder as a function of t. 
given the initial data 
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With the specified data we obtain the coefficients 

3 



1 2 , m 
%2 

,2,m 



(t) = 16 - 19* + 12t 2 - 3t 3 + -(1 - £) 4 atanh(i), 



(38) 



J l,2 

i2.m 
%2 
i2.m 
%2 

i2,m 
%2 



2(4 + t- 6t 2 + 3t 3 + 3(1 - <) 3 (1 + t) atanh(t)), 
(t) = \/6(5f - 3i 3 + 3(1 - tf{\ + tf atanh(t)), 
(t) = 2(-4 + f + 6i 2 + 3t 3 + 3(1 - t)(l + t) 3 atanh(t)), 

(t) = -16- 19t — 12t 2 - 3t 3 + 3 (1 + i) 4 atanh(t). 



In order to check how the code can cope with the singular behaviour at r = we 
use the generated numerical solution in the computational domain with < t < 1. 
In view of the expansion of the field components (f>k near r = we compute at every 
time-step 



(39) 



1 



0fc':T(*) = o#"-<M*,0). 



We evaluate the second derivative by using at least 4th order accurate one-sided 
finite difference operators which can be found in |37| and compare the result with 
the exact solution given in (38). Fig. [6] shows the numerically computed coefficient 




FIGURE 6. The first singular coefficient (f>±™ along r — for times 
0<t<l. 



04'™ and its exact value over the entire time interval [0, 1]. In Fig.|7]we zoom in to 
the last stages of the evolution and present the relative error for times larger than 
0.9. We see that we get roughly the expected 4th order convergence even though 
the accuracy is not great. This is clearly due to the numerical differentiation in the 
radial direction to obtain the coefficient in front of r 2 . 
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Relative error on the cylinder 
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Figure 7. Relative error in the computation of the first singular 
coefficient ^4'™ along r — for times 0.98 < t < 1. 



5. Discussion 

In this paper we have presented tests and preliminary applications for the lin- 
earised general conformal field equations in a setting near space-like infinity i of 
Minkowski space-time. We have used Friedrich's conformal representation of a 
neighbourhood of i° which separates the 'end points' of J^± by a cylinder /. 

Our motivation for this is the possibility to study entire space-times on a finite 
grid without the introduction of artificial boundaries. In contrast to traditional 
numerical approaches with a boundary at a finite distance from the source we are 
guaranteed that the corresponding solution is indeed asymptotically flat. In this 
setting we have access to asymptotic quantities such as the ADM mass and the 
radiation field at infinity. Furthermore, it is suggested by Friedrich's analytical 
results that, in principle, we have complete control over the smoothness of null 
infinity by choosing the behaviour of initial data on 1°. This issue is of fundamental 
importance, since gravitational radiation is defined unambiguously if null infinity 
has a certain level of smoothness. However, there still exist several outstanding 
questions even on the analytical side, and we hope that ultimately our numerical 
work may shed further light on these. 

The purpose of the present paper is to study the behaviour of the general con- 
formal field equations near the cylinder representing i° and, in particular, near the 
critical sets where touches the cylinder /. We limit our discussion to two 
particular conformal representations: one in which J> 'is diagonal', i.e., it is always 
transverse to the t = const, space-like hyper-surfaces; and the 'horizontal repres- 
entation', in which the t — const, hyper-surfaces tilt and become null for t = ±1, 
coinciding with J? ± . In an (r, t) diagram this is reflected by horizontal J? ± , as 
shown in Fig. [T] 
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Our numerical system was chosen as simple as possible. Apart from the linearisa- 
tion we also made use of the spherical symmetry of Minkowski space to decompose 
the spin-2 field representing the perturbed (rescaled) Weyl curvature into harmonic 
modes. 

Using this system we have been able to demonstrate that in the diagonal rep- 
resentation we can reach the singular set / + at t = 1 without loss of convergence 
in several test cases. This is not so surprising since the equations become singular 
'only' at /+, i.e., (r, t) — (0,1). With the explicit RK4 method that we use we 
never actually evaluate the equation at t = 1 so we never see the singular coeffi- 
cient 1 — k'(0) in front of dt4>i- Hence, the loss of hyperbolicity at I + does not 
seem to be important. However, as mentioned already in sect. |4.2| the break-down 
of hyperbolicity beyond t — 1 manifests itself very quickly. 

The situation is different in the horizontal representation. Here we cannot reach 
t = 1 because the equation for 04 becomes singular for all values of r and the 
characteristic speed of 04 blows up at t = 1. This has the consequence that running 
the code with any fixed time-step inadvertently leads to instabilities near t = 1. 
We can get around this by using adaptive time-stepping with the consequence that 
we can get arbitrarily close to t = 1 but in arbitrarily long wall time. 

There are several questions which need to and will be addressed in the future: 
In the diagonal representation we cannot compute far beyond the singular set. As 



we pointed out in sect. 4.2 the reason is due to the loss of hyperbolicity in every 
neighbourhood of I + . Consequently, the code is very unforgiving and blows up very 
quickly. If we were able to successfully compute the numerical solutions beyond the 
critical sets, this would mean that we could choose a hyperboloidal surface and 
obtain the induced data from the numerical solution. It would then be possible 
to use those as initial data for other numerical hyperboloidal codes, for example 
[18j [HJ [15J [HI Ell EH EH]. This would allow us to compute a large physically 
relevant portion of null infinity - and therefore the gravitational wave forms - in 
an unambiguous way. 

In future work we intend to explore the various possibilities of numerically con- 
structing hyperboloidal hyper-surfaces. Apart from the attempts to compute bey- 
ond the critical sets into the unphysical manifold we have started to implement the 
equations without the mode decomposition and will report on this work elsewhere. 

Another aspect that needs to be explored is the use of a compact conformal 
representation such as the one which exhibits Minkowski space as a subset of the 
Einstein cylinder. This would allow us to eliminate the artificial boundary at r = 1. 
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Appendix A. An exact solution of the spin-2 equation 



We very briefly describe how to derive an exact solution of the system ( 15 1-( 16 1. 
A detailed derivation can be found in |40j . The idea is to produce a solution on 
Minkowski space and then trace through the transformations described in sect. |2.1 
to obtain a solution in the conformal representation where spatial infinity is rep- 
resented as a cylinder. 
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The spin-2 system on Minkowski space with the metric g ([I]) using modal de- 
composition consists of the evolution system 

Rd T ^o - Rd R (j)o - O + a 2 0i = 0, 
2Rd T (t>i + 20\ - a 2 o + a 4> 2 = 0, 

(40) 2i?9 T 2 -a o (0i-0 3 )=O, 

2i?<9r03 — 203 — «O02 + C*204 = 0, 

RdT<f>4 + Rdn4>i + 4 - a 2 03 = 0, 

and the constraints 

2Rd R (f>i + 60i - a o 2 - a 2 0o = 0, 

(41) 2Rd R 4> 2 + 602 - ao03 - "001 = 0, 

2Rdn<f> 3 + 60 3 - a o 2 - «204 = 0. 

The first two evolution equations and the first constraint equation can be used to 
obtain a 2nd order wave equation for 0o alone: 

(42) i? 2 a|0 o -i? 2 ^0o + 4i?a T 0o-6i?9K0o + (a2-4)0o = 0. 
With the ansatz 

fa(T,R) = K{R)T(T) 



eq. (42) takes the 'almost' separable form 



T T K" K' 
R^+A^-=R—+Q 



T T K 11 R ' 

where / and /' denote differentiation of a function / with respect to the time and 
spatial coordinate, respectively. On the right hand side we have a function of R only 
while the left hand side is a linear polynomial in R with T-dependent coefficients. 
Taking another i?-derivative of this equation yields an equation between a function 
of R on the right and a function of T on the left. Thus, they must be constant. 
This leads to the ansatz: 

7?" TZ' a 2 - 4 

(43) fl | +6 |"V = W + 4m ' 

where the quantities appearing on the r.h.s must be arbitrary complex constants. 
Alternatively, the above ansatz can be interpreted as 

T T 

(44) R^-+A^-=kR + Am. 



Obviously, one can set 



T T 



T -k T 

r ' r 



Thus, the most general solution of (44 1 is of the form 
(45) T(T) = e mT , 

up to an irrelevant scale factor. The equation for the spatial part also admits an 
analytic solution. In order to keep things as simply as possible, we will consider 
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2:5 



here only the first non-trivial case 1 — 2. For this choice of Z (43) admits a solution 
of the form 



1 



(46) 7Z(R) = (e mR (3 - 6mR + 6m 2 R 2 - 4m 3 i? 3 + 2m 4 i? 4 ) c x - 3e" mB c 2 ) 



Using (45 I, (46) and the first four evolution equations, it is straightforward to show 
that 
(47) 





f e m{T+R) 








' e va{T+R) 








Lm(T+R) 


R 5 






' e va{T+R) 








f e m{T+R) 







IC 2 



}_ ^ e m(T+H)(_ 2 + mR) ^ + e m(T-fl)^g + g mR + 6m 2 R 2 + 2 m^R i ) C 2 ) 

d + e m(T_ii) (3 + 6mi? + 6m 2 i? 2 + 4m 3 i? 3 + 2m 4 i? 4 ) c 2 



Relations (47 1 provide a family of solutions for the system (40 1— (41 ), which satisfy 
a separation of variables ansatz. This family describes a superposition of spherical 
ingoing and outgoing waves. Before we continue to the derivation of the solutions 
for the complete system (15 1— (16), we will make a specific choice for the constants 
ci,C2,m, namely c\ = \,c% = 0, m = 0. With this choice (47) reduces to a static 
solution: 



(48) 



»0 = 



1 

#5' 



2 



»2 = 



V6 
i? 5 ' 



2 



1 



Using this static solution we now perform successively the transformations de- 
scribed above to obtain the corresponding solution with respect to the metric 
g, see eq. Q5J) . The main ingredient we have to be aware of is that according 
to [27] the spin-2 zero-rest-mass field transforms, under rescalings of the form Q, 
as 4>abcd = 0" 1 4>abcd- Using this formula and the appropriate rescalings of the 
spinor components and choosing the conformal representation with n(r) = r/(l + r) 



we, finally, end up with an exact solution of the complete system (|15|-(|16|: 

<l>o(t,r] 



(49) 



<p2(t,r 
(f>4,(t,r 



r 2 (l + r-t) 4 
(1 + r) 7 ' 
2r 2 (l + r - tf(l + r + t) 
{l + rf ' 
\/6r 2 (l + r-<) 2 (l + r + <) 2 
(TT7)7 : 
2r 2 (l + r-t)(l + r + <) 3 

r 2 (l + r + f) 4 
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Appendix B. Spin-weighted spherical harmonics 

We briefly summarize the derivation and conventions of the spin-weighted spher- 
ical harmonics s Y/ m from |27| and compare them to Friedrich's definition of the 
functions Ti 3 k , see for example |25j. 

Let V 1 be the (I + l)-dimensional vector space of complex homogeneous polyno- 
mials of degree I in two complex variables z\ and Z2, i.e. the vector space spanned 
by the basis {(f> l , . . . , (f>\} where 4>\,(z\,Z2) := z^z^ ■ Let a sesquilinear form be 
given on V 1 by 

(<&,4? m ) :=k\(l-k)\6 km , 

where our convention is that this form is anti-linear in the second argument; this is 
a scalar product on V 1 . For every t £ SU(2) in the standard matrix representation 



(50) t = 



i o 1 1 

t Q t 1 



with complex numbers t° Q , t i, t l $ and i 1 ! satisfying 



t ot i — t it o - i, t o-t i, t i 



we consider the map 

T t :C 2 ^C 2 , (zo^if^i' (^o^i) T , 

where "•" denotes a matrix product. For every non-negative integer I, this induces 
a representation Er of SU (2) on as follows 

U l : SU(2) x V 1 -> (t,</») ^ [/ t V := 4>oT t -i. 

It is shown in |41j that each of these representations is unitary (up to normalizations 
with respect to the above scalar product, see below) and irreducible. Indeed every 
irreducible unitary representation of SU (2) is equivalent to U l for some choice of I. 
For every 1 = 0,1,..., and to, k = 0, . . . , I, and t € SU(2), we find 

uW h ) = m\{l - m)\ (J) t^ [bl t a \ 2 ■ ■ ■ t°'K l)fc) 



771 ' 



where the notation (a\a2 ■ ■ ■ a{) m means the following: (i) symmetrize over the 
indices a\, . . . , aj, (ii) replace to indices by and I — to indices by 1. A normaliza- 
tion yields the matrix elements of the unitary representations of SU (2) ; those are 
Friedrich's functions 

^kit) ,,„„,' „ 777^7 ' ^ f ''<•'■■ 

(51) 

According to the Peter- Weyl theorem |U] the functions \Jl + 1 Ti m k , with I = 
0, 1, 2, . . ., and m,k = 0,1, ... ,1 form an orthonormal basis of the Hilbert space 
L 2 (SU(2)) with respect to the standard normalized Haar measure. Hence, every 
function on SU(2) can be expanded in this basis. 
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Now, for any non-negative integeijjZ, let the space W 21 be the (2/+l)-dimensional 
vector space of those spinorial tensor fields on Minkowski space which are spanned 
by the basis elements 

Z(l, m)A 1 ...A t - m B 1 ...Bi +m '■— °(Ai ' ' ' OAi_ m ^Bi ' ' ' is l+m )j 

for m = — I, . . . , I, see Eq. (4.15.93) in [27]. Here (b A , i A ) is the "constant dyad" of 
Penrose (as opposed to a "rotated dyad" below) which is normalized by 1 = oaX ■ 
Clearly there is an isomorphism between '5 : W 21 — > V 21 with 

* : Z(l, m) Al ... A2l ^ Z(l, m) Al ... A2l Tl A i . . . IP 4 - = cf>fL m , 

for 

II A := -z x b A + z l A . 

Roughly speaking, each oa is replaced by zq and each I a by z\. By using this 
isomorphism, we can induce a scalar product on W 21 from the one above on V 21 
(in the following we do often not write spinor indices when it is convenient): 

(Z(l, m), Z(l, k)) := m)), fc))> = (/ - fc)!(Z + fc)!<W 

However, we also have a geometric scalar product on this space of space spinors in 
Minkowski space constructed as follows. Let 

T AA ' := l = (d A o A ' +t A l A '). 

Phen we set 

(Z(l, m), Z(l, k)) G := T AlA>1 ■ ■ ■ T A ^Z(l, m) Al ...A ai W^)Ai L ...A> 21 ■ 

Noting that T AA 'b A = l A ' /V? and T AA ' l A = -b A ' /y/2, we find 

{l -k)\{l + k)\ _ 1 
2' (2/)! mfc ~ 2 ; (20! 

cf. Eq. (4.15.94) in [27]. 

We can also use the isomorphism \I> and the representations U 21 of SU(2) on V" 2 ' 
above to define representations of SU{2) on W 21 , which we also call U 21 , as follows: 

U 21 : SU{2) x W 21 -> TP 2 ', (t,Z(l,m)) h> U 2l Z(l,m) := ^~ l {Uf V(Z(l, m))). 

We can write this as 

Ut Z(l,m)A 1 ...A l _ m B 1 ...B, +m = ' ■■0 Al _ m LB 1 • ' ' ^B !+m ) 5 

where 

(52) oa = t° 6,4 + t^M, ^ = t°iOA + 

In (4.15.95) of [27], one defines the functions 

,Zj, m (t) := m)^...^.^..^^ • • • t ^- ^-+i . . . o Ml \ 



(Z(l, m), Z(l, fc)) G = V /; V 7 ^ = (Z(l, m), fc)> 



'^All of the following constructions also work for half-integer values of I. In this case, the other 
indices m and s used in the following must also be half-integer valued. In this paper, we focus 
exclusively on the integer case. 
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on SU(2) (the quantity W in (4.15.95) in [27] is identically unity in our case of 
pure rotations). Noting that T AA ' oa = l a ' /V% and T AA ' la = -o A ' /\f2, we can 
write this as 

s z hm (t) = (-ly^^T^^.-.r^^za^)^...^^^/^)^..^ 

= (-l) l +°2 l (Z(l, m), Uf l Z(l, s)) G = t^L ( Z (l, m), U?Z{1, a)) 

It follows 

s Zi, m (t) = i -^^ T ^ m is( t ) y/(l-m)Kl + m)l(l- S y.(l + S )l 

The spin-weighted spherical harmonics s Y/ m are defined by (4.15.98) in [27] as 
normalizations of the functions s Z^ m (t): 



.rut) -(-ir\W 



4tt(Z - m)!(Z + ra)!(Z - s)!(Z + s)f 



Hence, it follows 



(53) s Y iro (t) = (-l)*+™^±! T 2l l - m l _ s {t). 

We see explicitly that these functions are defined for every I = 0, 1, . . . and for every 
m and s constrained to m — —I, . . . , I and s = —I, . . . , I. 

The reader should notice that, while the representation spaces V for Z = 0, 1, . . . 
give rise to aZZ unitary irreducible representation of SU (2) and hence to a basis 
of the HUbcrt space L 2 (SU(2)) as above, the representation spaces IT" 2 ' (which 
are isomorphic to V 21 ) for integer- valued I give rise to only "half" of the irreducible 
representations. In particular, it follows that the functions s Yi m with integer-valued 
indices Z, m and s do not form a basis of L 2 (SU(2)). However, for a fixed value 
of s. the functions s Yj m form a basis of functions with spin-weight s in L 2 (SU(2)). 
This is not only true for the integer spin case, which we have considered here, but 
also for the half-integer spin case. 

So far, we have considered s Yim as functions on SU(2), and not, as it is usually 
done, as functions on S 2 . We argue as follows without going into the details (which 
can be found for instance in [42 ). It is a fact that SU{2) is diffeomorphic to S 3 , and 
S 3 is the Hopf bundle over S 2 with the Hopf map 7r as the bundle projection. This 
bundle is not trivial and hence smooth functions on SU (2) do in general not give 
rise to smooth functions on S 2 . We can, however, choose a smooth local section 
in the bundle, i.e. a smooth map a : U — > SU(2) compatible with n where U is 
an open subset of § 2 . With this the smooth functions s Yi m can be pulled back 
to smooth functions s Yi m o a defined on U. The choice of local section made by 
Penrose is the following. We introduce standard polar coordinates (9, <fi) on S 2 and 
set U :— S 2 \{6* = 0,7r}. Then, the following defines a smooth local section: 

tt OTTf^ to ^ * f t0 o tQ i\ /V^cosf -e-^/ 2 sin§\ 
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The reader should compare this to Eq. (4.15.123) in [27] 

t i t o 
i 





b A o A o A i A \ _ ( e-^/ 2 sin | a" 1 */ 2 cos § 
£ A o A J ~ l -e^/ 2 cos I e^/ 2 sin e 



2 ° kJiii 2 

which is consistent. Using these coordinate expressions, we can compute the func- 
tions Ti m k as functions on S 2 via Eq. (511 and then s Y/ m by Eq. (53 1. It is im- 
portant to notice that while T™ k and s Yj m are smooth functions on SU(2) and 
their pull-backs to U are smooth, there is in general no smooth (not even a continu- 
ous) extension from the dense subset U to S 2 . Here are some examples of explicit 
expressions for some harmonics 



0^1,-1 — 



sm ae 



o*Lo — 



COSt 



0*1,1 — 



sm( 



o ir t> 



and 



>r 2 



, sin 4 ^e- 21 
4tt 2 



2^2.0 



2*2,-1 



, sin 2 - sin 9e 
4tt 2 



sm 0, 

32tt 



etc. 
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